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1 . STATEMENT O F WORK 


The statement of work includes data analysis and supporting research 
in connection with the following broad objectives: 

(1) Provide a precise and accurate geometric description of the 
earth’s surface. 

(2) Provide a precise and accurate mathematical description of 
the earth's gravitational field. 

(3) Determine time variations of the geometry of the ocean sur- 
face, the solid earth, the gravity field and other geophysical 
parameters. 
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2. ACTIVITIES RELATED TO THE NGSP 
(Grant No. NGL 36-008-093) 

2.1 Data Acquisition and Processing 

The data of the WEST (West European Satellite Triangulation) and the ISAGEX 
(International Satellite Geodesy Experiment) programs are at our disposal . The pur- 
pose of this investigation is to utilize some or all of the above observations in order 
to improve the values of some station coordinates on the European continent which 
are presently included in the WN-14 solution and to assess the quality of the WN-14 
solution with the help of the new data available. A detailed description of the data 
is given in the previous Semiannual Status Report. The current status of acquisi- 
tion and processing is given below. 

2.11 WEST Data 

There are two sets of optical data available . One set contains the direction 
cosines of single ficticious images per plate including the standard deviations which 
were derived from polynomial fitting. The other set contains the direction cosines of 
seven fictitious images. All directions are given in the Greenwich Hour Angle 
Declination system. 

2.11.1 Single Image Data Processing 

Since no program was available at OSU to process single image data, the 
,OSUGOP program (which was previously used to process BC-4 seven image data) 
was suitably modified by James P. Reilly. The subroutines READIN and ASD 360 
had to 1-3 completely rewritten, while the subroutine FORMRN had to undergo only 
minor changes . . 

Transformation of variances : The variances of the observations were given in 

the form of standard deviations along and across the satellite trail. The modified sub- 
routines require as input the standard deviation of the Greenwich hour angle multiplied 
with the cosine of the declination, the standard deviation of the declination and the 
covariance term. The variances were transformed as follows: 
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Fig. 1 Transformation of Variances 


AL denotes length of the trail 

6 a » denote declination of satellite (beginning and end of trail) 

Ot denotes rotation angle . 

The actual rotation is approximated by a rotation around the point A where the 
satellite trail AB is taken to be a straight line. 



Fig . 2 Rotation to GHA System 


We obtain the relation 


/x\ / cos(-g) sin(-O') wx'v 

\y/ y-sin(-a) cos(-a) Ay/ 


cos(-o;) sin(-O') 
a) cos(- 

where a is computed from the spherical relation 


sin ot = 


sin(6 E - 6 a ) 
sin AL 


( 1 ) 


( 2 ) 
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Using the given variance-covariance matrix 

s.v = ( 1 ( 

\o 6 y ?/ 

where 6* is the variance along the trail and 5 2 is the variance across the 
trail, we obtain with relation (1) and the law for propagating the variances, the 
transformed variance-covariance matrix 



'cos" (-a)5 2 + sin 2 ( - a )6 c', sin(-o:)cos(-Q!)(-6 2 + 6 2 )) 


' sin (-a)6 2 + cos 2 (-a)6| 


It should be noted that 6 x denotes 6 qHA • cos 6 and that 6 y denotes 6 6 . 

Some information about the observations stations : Station numbering : 

The whole set of observations contained 30 different tracking stations. They are 
listed in Table 1 and their relative location can be seen from Figure 3. Since the 
numbering system for the WEST stations and WN-14 stations are independent, the 
Same station number was assigned to different stations. In order to avoid confusion 
the WEST stations were in part renumbered with a four digit number . The complete 
station number consists of four digits, where the first two digits were arbitrarily 
chosen as 87 and the last two digits were the same as in the corresponding WEST station 
number. In some cases even the second digit had to be changed. The following 
table presents both numbers . Further modification of the station numbers might 
become necessary . 
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Table 1 Summary of Observation Stations, WEST 


Relative constraints for nearby stations : For the adjustment computations it 
is important to establish the exact relationship between nearby stations which can be 
introduced as relative constraints. In most cases the following information was 
extracted from Circular Letters which were distributed by the International Assoc- 
iation of Geodesy during the time of the WEST campaign. In some cases the relative 
location of observation stations could be established from Cartesian coordinates 
which are given in Ehmsperger [4] . The coordinates are given up to millimeter 
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and it seems reasonable to assume that they accurately render the relative position 
of observation stations. (The actual distance between nearby stations usually amounts 
to a few meters.) The details about some specific stations and our observations and 
conclusions are given below. 

STATIONS 8702 (10002) MADRID - 8703 (10003) MAORI 
Circular Letter No. 35: The new pillar 10003 MAORI (IGN camera) is 2.74m from 
10002 MADRD. Ehmsperger [4] gives the following coordinate differences: Ax = 

-0.065, Ay 2 . 696, Az 0 , 282. Check: the linear distance as computed from the 
coordinate differences is 2.71m. Conclusion: the above coordinate differences 
seem to be correct. 

STATIONS 9001 DELFT - 9002 DELFY 

Ehmsperger [4] gives the following coordinate differences: Ax = 3709. 057, Ay = 
1053.539, Az = 2925. 820. Both stations are listed at the NASA Directory of Observation 
Stati. n Locations. Their differences agree with the above values. 

STATIONS 6010 HOPBG - 6110 HOPBG 

Circular Letter No. 31 confirms that both stations are identical. The change of 
the station numbers agrees with the convention during the WEST campaign that a 
change of the camera should be indicated by a change in the third digit of the station 
number. In this case station 6010 HOPBG was equipped with an IGN camera while 
6110 HOPBG carried a BC-4 camera. 

STATIONS 6004 BRNSG - 6012 WSNDF 

Circular Letter No. ^3: both stations are approximately 2m apart. Ehmsperger 
(4] gives the following coordinate differences: Ax = 1.785, Ay = - 1. 335, Az = - 0. 161. 
Check: the linear distance as computed from the coordinate differences is 2. 23 m. 

The coordinates of station 6004 BRNSG are listed in the NASA Directory of Observa- 
tion Station Locations and agree with those used by Ehmsperger. Conclusion: the above 
coordinate differences seem to be correct. 

STATIONS 8008 TANIA - 8709 CATNA 

Ehmsperger [4] gives the following coordinate differences: Ax = 13333.95, Ay = 
10080.89, Az = -22966.80. 

STATIONS 8004 CATAN - 8008 TANIA 

Circular Letter No. 37: the new station 8008 TANIA is approximately 12m south, 
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southesat of station 8004 CATAN. Ehrnsperger gives the following coordinate differ- 
ences: Ax= -4.036, Ay = -8.238. Az - 7.882. Check: the distance as computed 
from the coordinate differences is 12.06m and points in the expected direction. 
Conclusion: the above coordinate differences seem to be correct. 

OLD BC-4 SITE 6016 CATANIA 

Circlular Letter No. 37: the old BC-4 site is 2.76m south of station 8004 CATAN. 

This information made the computation of relative constraints possible (assuming that 
both stations have the same heights): Ax = -1.67, Ay = 0.00, Az = 2.17. 

OLD BC-4 SITE 6065 HOHENPEISSENBERG 
WEST observation station 6010 HOPBG is identical with station 8032 MUNICH 
according to [6]. The following relative constraint between the WN-14 station, 6065 
HOHENPEISSENBERG, and the WEST station, 6010 HOPBG, were derived from 
geodetic coordinates given in [6]: Ax= 21.26, Ay = 54.46, Az = -24.52. 

In order to compare or to combine the WN-14 and WEST systems, common 
stations have to be identified. The following table of identical stations could be gathered: 


WN-14 No. 

Name 

WEST No. 

Reference 

i 6006 • 

TROMSO 

8742 

[7] No. 26 

6016 

CATANIA 

8004* 

[7] No. 37 

6065 

HOPBG 

6110, 6010* 

[6] 

•8009 

DELFT 

9001 

[6] 

8010 

ZMDLD 

8721 

[6] 

8011 

MLVRN 

8732 

m 

8019 \ 

NICE 

5004 

[6] 

8030 

MEUDN 

5001 

[6] 


* 

See the specific information given for the old BC-4 sites 
Table 2 Identical Stations for WN-14 and West 


Preliminary Adjustment. Computations 
Adjustment WEST No . 1 

The purpose of the first adjustment was to find the adjusted variance of unit 
weight and to get a first insight into the quality of the data. The following input data 
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and constraints were considered: 

1) The transformed variance-covariance matrix was used as described 
previously. 

2) For all identical stations as given in Table 2, the approximate coor- 
dinates of the WN-14 solution were used. All other station coordinates were trans- 
formed from the European Datum EU-50 to the OSU (WN-14) datum using the following 
parameters: 

EU-50: A = 6378388.0m, 1/F = 297.0 

OSU (WN-14): A = 6378155.0m, 1/F - 298.2494985 

dx = -99.4, dy = -132. 0, da = -116.0, ds = 0 .0675 ppm. 

3) AH relative constraints which were previously described were 
enforced by appropriate weighting. The weights are based on an assumed accuracy 
for the geodetic survey of approximately 1:50,000. 

4) The origin of the cooidinate system was defined by Inner Adjustment. 

5) The scale was introduced through the base line 6016 CATANIA - 
8742 TROMSO with an accuracy of 1 ppm. 

Result of the adjustment: 

a) A posteriori variance of unit weight: 35.5. 

b) In all cases the adjusted coordinate differences of those stations which 
were connected by relative constraints do not deviate from the the constrained values. 

Adjustment WEST No. 2 

The purpose of this adjustment is to combine the WN-14 system and the WEST 
system. The adjustment is based on the following input: 

1) All possible relative constraints. 

2) Chord constraints between stations 6016 CATANIA and 8742 TROMSO. 

3) The coordinates of all common stations between both systems 
(Table 2) are constrained to the adjusted WN-14 values using weights which were 
computed from corresponding variances. 

4) The heights of all common stations are constrained to the values 
given in [3], Table 3. 3-3. 

Result of the adjustment: a) A posteriori variance of unit weight: 37. 1. 
b) Comparison of coordinates of identical stations: 
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Station No. 

Ax 

(m) 

CTx WN-14 

1 

ay WN-14 

m 

a 2 WN- 14 

WN-14 

WEST 

Ox WEST 

Oy WEST 

__ 

°z WEST 

6065 


-0.56 

2.02 

1.83 

HU 

mmzm 

0. 81 

wbbbm 

6016 


■ 

1.81 

1.67 

1.47 

mm 

0.07 

2.24 

2.08 

6006 

8742 

-0.35 

2.36 

2.16 

3.03 

2.92 

2.58 

-0.94 

2.89 

2.31 

8030 

5001 

■ 

6.46 

4.94 


mm 

BB 

0.79 

1. " 

8019 

5004 

0.69 

4.12 

3.50 


7.91 

6.64 

0.38 

4.31 

3.68 

8010 

8721 

■l 


■ 

8.28 

5.02 

-4.. 89 

5.44 

3.57 

8011 

8732 

m 

8.86 

3.61 

-36.81 

14.27 

5.31 

4.48 

6.96 

3.84 

8009 

9001 

■ 

8.48 

4.28 

-6.92 

Kfe $ 

-0.07 

6.8b 

4.10 


Table 3 Adjusted WN-14 Coordinates Minus the 
Coordinates of Adjusted West No. 2 

c) In Tables 4 and 5 the results are given for a transformation of the 

coordinates obtained from adjustment WEST 2 to the WN-14 system. 

2.11.2 Seven Image Data Processing 

The seven image data received on cards was first transferred onto tape. 

This data gives direction cosines, event/stationwise. Some event numbers have 
been duplicated and there are some image numbers from 8 to 14. A second program 
read the data from the tape, transformed the direction cosines to right ascension 
and declination, generated additional parameters (such as numbers of stations in each 
event) required for input for the OSUGOP program and transferred the data in the 
new form onto disk. A third program has been made to read the single image data 
(for variances) from cards, match it with the modified seven image data (event and 
stationwise) and output the merged observational data giving a, 6 and variances for 
the seven images. A variational parameter has been added to vary the variances 
of images 1, 2, 3,5, 6, 7 with respect to the variance of the 4th image which is available 
in the single image data . 
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SOtUriON“FOR~3 'TRANSrATION~l“SCAC!=" AND “3" ROTATION "PARAMETERS 

COSING VARIANCES ONLY) 

DX Dy Dz DELTA OMEGA PSI EPSILON 

METERS METERS METERS (XI. 0+6) SECONDS SECONDS SECONDS 

"-TV24 -1~T2? -OT6D 0^51 0^29 -0V01 -0 .18 

+ 2.57 + 2.12 ± 3.12 ± 1.05 ± 0.52 ± 0.26 ♦ 0.67 


VARIANCE -"COVARIANCE MATRIX 


S02 = 

0.79 





0.659D+01 

0.159D+01 

0.493D+01 -0.241D-07 

-0.558D-05 

0.830D-06 

0. 5780-05 

0.1590+01 

0.4510+01 

0. 2240+01 -0.1090-05 

—0. 185D-05 

0.6600-06 

0.296D-05 

0.4930+01 

0.2240+01 

0.975D+01 0.219D— 06 

-0. 576D-05 

0.142D— 05 

0.8950-05 

-0. 2410-07 

-0.109D-05 

0. 2.19D— 06 0.11 ID— 11 

0.179D-13 

-0.1850-12 

-0. 8800-13 

-0. 5580-05 

— 0.185D-05 

-0.576D-05 0.179D-13 

0.643D— 11 

-0.1400-11 

-0.6750-1 1 

0. 8300-06 

0. 6600-06 

0. 1420-05 -0.1850-12 

—0. 140D-1 1 

0.159D-11 

0.177D-11 

0.5780-05 

0.296D— 05 

0.8950-05 -0.880D— 13 

* 

^ Jr 

— 0.675D— 11 

0.177D-1 1 

0. 1050-10 



m+HlltNIb Ur LUKK tLA 1 1 UN 



0. 100D+01 . 

0.2920+00 

0.6160+00 -0.893D-02 

-0.8580+00 

0.2560+00 

0.6950+00 

0.2920+00 

O.IOOD+Ol 

0.338D+00 -0.4800+00 

-0.3430+00 

0.2460+00 

0.4310+00 

0.616D+00 

0.338D+00 

0 • 100D+01 0.668D-01 

-0.7270+00 

0.360D+00 

0 • 884D+00 

—0.8930-02 

— 0.488D+00 

0.668D-01 0 .1000+01 

0.6710-02 

—0. 139D+00 

-0.258D-01 

— 0.858D+00 

— 0.343D+00 

-0.727D+00 0*67 ID— 02 

0.1000+01 

-0.4370+00 

-0.8220+00 

0.256D+00 

0.2460+00 

0 • 360D + 00 -0.139D+00 

— 0.437D+00 

0. 100D+01 

0.432D+00 

0.6950+00 

0.431D+00 

0. 884D+00 -0.258D-01 

-0.822D+00 

0.432D+00 

0. 1000+01 


Table 4 


RTSTOTJATS 


DX DY D2T 


600 6 - 

_ T _ T .- 

-TT9 

? ^ 1 

6016 

1.8 

-2.1 

-O.l 

606 5 

' . Am*, h 

-1.3 

-0.8 

8009 


5.5 

-0.4 

80110 



4.5 

8011 


EL • il 

-5.6 

8019 


Bli ■ ;■ 

-1.0 

8030 

1 bv % j 

Bf. ■ H 

-1.7 


Table 5 
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After outputting the merged data onto a tape in a form compatible with the 
input requirement of OSUGOP, experiments similar to the ones done with the single 
image data will be performed. 

2. 12 ISAGEX Data 
General Remarks 

A detailed description of the ISAGEX data as obtained from the Centre National 
d'Etudes Spatiales is given in the previous Semiannual Status Report on pages 4-10. 
The data consist of laser ranges and optical observations. 

It was already reported that no simultaneous laser range observations could 
be found. Therefore, efforts were made to further process the optical observations 
only. The preliminary results, which were already reported, indicated either a 
very poor quality of the ISAGEX data or a blunder in processing. It was therefore 
decided to completely re-examine the investigations done so far, starting with the 
preprocessed data as provided by Wolf Research and Development Corporation. 

The data did not include observations of MIDAS 4 and PAGEOS. These data were 
thus sent to Wc^f Research and Development Corporation for preprocessing. However 
it was learned during a telephone conversation with Ms. Donna Walls of the Wolf 
Corporation, that difficulties had arisen in obtaining the correct input data for pre- 
processing. It was consequently decided not to use the MIDAS 4 and PAGEOS data 
anymore for this investigation. 

Current Status of Processing the ISAGEX Optical Data 
The preprocessed data were tested for simultaneity. Allowing a time gap of 0, 2ms, 
a total of 353 observations proved to be simultaneous, involving 13 different stations 
which are exclusively located in Europe. 

hi the next step the quality of the observations was tested by forming the 
normal equations using the OSUGOP program. As explained in the Reports of the 
Department of Geodetic Science No. 190, page 12, a so-called test distance can be 
used to specify rejection criteria for each observation; or, conversely, by 
looking at the computed test distance we can judge the quality of the observation. 

A large test distance indicates either bad quality of observations or poor approxi- 
mate station coordinates. In our computations we used first, station coordinates 
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in the European system ED 50 as extracted from the NASA Directory of Observation 
Station Locations and from various ISAGEX documents, and secondly, the approximate 
coordinates which were used for the WN-14 solution. If no WN-14 coordinates were 
available, approximately transformed coordinates were used. The coordinates of 
station 1147 ONDREJOV were extracted from [5]. Both computations showed only 
minor differences. In Figure 4 the distribution of the 13 remaining observation 
stations is given. Also the number of observations on each line with a test distance 
smaller than 15 arc sec are shown. It is important to note that only a few observa- 
tions which qualify, exhibit such a large test distance, while the vast majority has 
a test distance of 1 arc sec or even less. The remaining observations have in 
most cases, test distances of several thousands arc sec, which probably indicates 
an error in data reduction. A closer inspection of those lines shows that the 
approximate coordinates do not cause such large test distances . 

The following table lists the 13 qualifying stations: 


ISAGEX No. 

Name 

1055 

UZHGOROD 

1184 

RIGA 

1147 

ONDREJOV 

1181 

POTSDAM 

8009 

HAUTE PROVENCE 

8010 

ZIMMERWALD 

8011 

MALVERN 

8019 

NICE 

8031 

EARLY POINT 

8034 

YPBURG 

1 9004 

SAN FERNANDO 

9030 

DYONISOS 

9120 

SAN VITO 


Table 6 ISAGEX Stations 


Some Remarks on the Observation Stations : 

The two data sets, WEST and ISAGEX, are independent sets, but the observations 
were made from common or nearby stations. It is, therefore, important to uniquely 
identify the observation stations. This investigation is still in progress and the 
following two tables are only of preliminary character: 
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ISAGEX No. 


Name 

mam 

9004 

SAN FERNANDO 


8015 

HAUTE PROVENCE 

8010 

8010 

ZIMMERWALD 

8019 

8019 

NICE 

8011 

8011 

MALVERN 


Table 7 Apparently Common Stations between ISAGEX and WN-14 


| ISAGEX No. 

Name 

WEST No. 

Name 

j 8031 

EARLY POINT 

mmm 

EDNBG 

| 8034 

YPBURG 


YPBURG 

8010 

ZIMMERWALD 

■ 

ZIMMERWALD 

8019 

NICE 

5004 

NICEM 

5005 

GOULT 

8009 

HAUTE PROVENCE 


Table 8 Apparently Common Stations between ISAGEX and WEST 


Preliminary Adjustment Computations : 

In order to further test the observations, various adjustment computations were 
carried out. From Figure 4 it can be seen that the following four stations do not 
form closed figures: 

1055 UZHGOROD 
1084 RIGA 
9030 DYONISOS 
9120 SAN VITO. 

These stations have been neglected in subsequent computations. 

Adjustment: ISAGEX 1 

The adjustment is based on the following information: 

1) Number of stations: 9 

2) Standard deviation of the observation: 1 arc sec 

3) The origin of the coordinate system is defined by Inner Adjustment 

4) All five stations given in Table 7 have been constrained to the WN-14 
coordinates where the weights were computed from the standard deviations given 
in the Reports of the Department of Geodetic Science No. 199, pages 118-145, as 
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Weight = 1/a 2 . 

Result: 

a) Variance of unit weight: 4.0 

b) Comparison of coordinates: WN-14 - ISAGEX. 


Station No. 


Ax(m) 


Ay(m) 

.... 

Az(m) 


WN-14 

ISAGEX 




8010 

8010 

26.21 

■ 

-17.64 

8. 28 
5. 36 

-15.1 

5.44 

5.40 

8011 

8011 

2.43 

wmm 

-60.88 


BUM 

6.95 

9.04 

8015 

8009 

-28.01 

4.19 

6.36 

-21.47 

■ 

25.31 

4. 38 
6. 87 

8019 

8019 

■1 

4.12 

6.31 

44. 15 

mmm 

deal 

-12.45 

WMMM 

Rii91 


Table 9 WN-14 - Adjustment ISAGEX 1 


2.13 Determination of Transformation Parameters and Network 
Distortions 

The Fourteenth and Fifteenth Semiannual Status Reports contain various 
tables of transformation parameters and figures which give indications of network 
distortions. Investigations have been made so far for the North American Geodetic 
D>atum NAD 27 and the Australian Geodetic Datum. Both reports also contain a 
detailed description of the procedure used in this investigation. During the present 
reporting period, computations have been carried out for the South American Datum 
SAD 69 using Doppler stations which were provided by the Defense Mapping Agency 
(Attachment 1). The results are given in Table 10 and Figures 5, 6, 7, 8, 9 and 10 
At this time the investigation regarding the transformation parameters and 
network distortions using the previously mentioned procedure can be considered 
complete. It is intended to publish all results and updated computer programs in 
a final report. 
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NWL9 

SAD69 

Ax(m) 

-77. 8 ± 2. ^ 

Ay(m) 

-12.4 ± 3.9 

Az (m) 

-49. 5 ± 2.6 

A(10 -6 ) 

- 0. 99 ± 0.55 

wO 

1. 18± 0.33 

*<"> 

- 0 . 90 ± 0.13 

e(") 

0. 16± 0.10 

a (") 

0. 33 ± 0.10 

£ <") 

- 0 . 48 ± 0 . 10 

V(") 

- 1.37 ±0.35 


Table 10 Datum Transformation Parameters 
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Latitude Differences After Transforming SAD69 to NWL9D (in Metres) 
7 Transformation Parameters (NWL9D-SAD69) 


MODEL 


Molodenskii Veis 


Au (m)=-77.8 ±2.7 
Av (m)=-!2.4 ±3.9 
Aw(m) = -49.5 ±2.6 
A(ppm)= -0.99 ±0.55 


w( M )= 1 .18 ±0.33 
^(")=-0.90±0.l3 
c(")« 0. 16±0. 10 


a(")= 0.33 ±0.10 
£C)=-0.48±0.I0 
17 (")=- 1.37 ±0.35 






Longitude Differences After Transforming SAD69 to NWL9D (in Metre**) 
7 Transformation Parameters (NWL9D-SAD69) 


'MODEL 


Au (m)=-77.8 ±2.7 
Av (m) =- 12.4 ±3.9 
Aw (mt = -49.5 ±2.6 
A ( ppm) = -0.99 ± 0.55 


Moiodenskii 
w(")= 1.1 8 ±0.33 
^(")=-0.90±O.I3 
eC)= 0 . 1 6 ± 0. 10 


30121 


30196 


Veis 

a(")= 0.33 ±0.10 
£( W )=-0.48±0.I0 
1.37 ±0.35 
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Height Differences After Transforming SAD69 to NWL9D (in Metres) 
7 Transformation Parameters (NWL9D-SAD69) 


Au (m) =-77.8 ±2.7 
Av (m) = - 12.4 ±3.9 
Aw (m) 3 -49.5 ±2.6 
A (ppm) 3 -0.99 ±0.55 
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3. ACTIVITIES RELATED TO EOPAP 
(Grant No. NGR 36-008-204) 

3.1 Sea Level Slopes Along the Continental Boundaries of the U.S.A. 

The previous Semiannual Status Report sets out the statement of the problem 
and the conclusions reached in comparing the results of leveling as done by geodesists 
and oceanographers. These computations have been refined by using Prey reduction 
instead of free air reduction for gravity values. These modified computations, (see 
Attachment 2), do not change the conclusions reached earlier. 

The subject was discussed in Washington, D. C. on June 16, 1975 at the 
meeting of the American Geophysical Union by the Subcommittee on the Discrepancy 
in the Geoids. The following points emerged from these discussions : 

(i) The oceanographers indicated that further research effort in identi- 
fying the cause for discrepancy could be concentrated on the region very near the 
ocean surface. 

(ii) The oceanographers agreed to supply further details about their 
method of work by giving detailed calculations at one of the stations . 

(iii) The study could be extended to the continental coasts of both North 
and South America for which data is understood to be available . 

. The computations previously carried out were modified to depict the mean 
sea level with respect to the origin of geodetic leveling. This has been shown in 
Figure 11. This may be studied in reference to the previous Semiannual Status 
Report, Figures 3.1-3 and 3. 1-5. 

The graphics give no additional information. However, the agreement 
about the direction of slope between oceanic and geodetic leveling is now less 
obvious. Further investigations could be attempted after more data is received. 
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3 . 2 Rotation of the Earth 


The equations of motion of a rigid body about its center of mass are governed 
by the well known Euler's equation. Under the assumption of absolute rigidity, the 
relative positions of all mass particles constituting such bodies are constant, so the 
external form as well as the moments of inertia are fixed and independent of time. 

It appears clear now that the assumption of rigidity for the earth is incorrect 
and a better modeling should be investigated. Because the mass distribution of the 
earth is subject to variations with time (i.e. , tidal deformation, crustal motions, 
etc.), producing changes in its inertia tensor, the rotational dynamics of the earth 
are better studied by the Lagrange- Liouville equations. 

The mathematical theory can be summarized in the following three equations, 
expressed in matrix notation by 

[L] = [I] («& +[ij{«5+ {h} + [tO][I] {0)} + [w]{h} (1) 

where: 

{L} s vector of external torques 

{co 3 =• rotation vector of the earth 

{h] = relative angular momentum vector 

[I] = earth's inertia tensor of the second order 

[to] = skew-symmetric matrix of the rotation vector. 

All of the above vectors' components are referred to an arbitrary earth fixed system. 
Significant simplifications will be introduced when we choose the principal moments of 
inertia axis as the reference system. As usual, the derivatives respect to time are 
represented by a dot. 

Equation (1) can be solved for the general case of a deformable earth by 
taking into consideration the variations of [I] and [h]. For example, it is 
possible to write 

[I] = . [I 0 ] + [A I] P + [ AI h + [AI] t + [AI] E + other effects (2) 

where: 

[1 0 ] = initial value of the earth’s tensor of inertia 

[A I] p = contribution to [I 0 ] due to crustal mass displacements (plate tectonics) 
[AI] r = contribution to [I 0 ] due to rotational deformation 
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[All T = contribution to [I 0 ] due to tidal deformation 

[Al] E = contribution to [I 0 ] due to large earthquake faulting. 

Similarly, 

{h} = {h 0 }+ fh} P + fh} R + {h} T + {h} E + other effects . (3) 

Each of the above contributions to [I 0 ] and {h 0 } is obtained after consideration 
of the particular adopted earth model. 

For example, one may express: 

[AI] P = £[Al] P (4) 

1=1 i 

where n is the number of tectonic plates constituting the earth crust. 

Likewise, 

[AIJt S [AI] T * (5) 

j=i J 

where j - 1, 2, • • • m refers to the moon, sun and planets producing tidal 
deformations . 

Thus, the solution {to } of the differential equations (1) will provide the 
changes in the earth rotation vector after consideration of the latest geophysical 
theories. 

The nature of these changes will depend on the hypothesis about the 
distribution of masses and its time variations. Clearly, secular changes in 
{to} , if any, will be produced by [AI] P . Periodic changes will be caused by 
CAI) t and sudden variations in fu>} will be due to the effect of [Al] E . 

The present intent is to answer the controversial question of secular drift 
of the "mean pole. M The tectonic plate model given by [Solomon and Sleep] is used 
in the investigation for the computation of [A I] P (see equation (4) ). 

Analytical expressions for the obtention of [Al] P ' due to differential motion 
of the plates have been developed as well as formulas for the computation of 
(:Al] Pl and [h] Pj . 

A computer program is being written in order to obtain the contribution 
to [I 0 ] of each independent plate [Al] P . This program integrates over an 
ellipsoid and assumes the Heiskanen theory of isostatic compensation for the 
crust. This hypothesis roughly agrees with the separation crust- mantle as 
postulated by seismologists (i.e. , Mohorovicic or M discontinuity). 
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Also a computer program is under way in order to obtain good initial values 
of the earth inertia tensor [I 0 ] . 

In the future, after the values of [I 0 ], [Al] P * [Al] P , {h} P and {h} p are known, 
the differential equations (1) will be solved numerically, thus answering the question 
of a possible secular shift of the pole due to crustal unrest. 

REFERENCES 

Solomon, S. C. andN. Y. Sleep, 1974. "Some Simple Physical Models for Absolute 
Plate Motions," Journal Geophysical Research , 79 , 2557-2567. 


3.3 Close Grid Geodynamic Satellite (CLOGEOS) System 

Some effort was made to help the investigations performed for Marshall 
Space Hight Center. The results are included in Attachments 3 and 4. 
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147 Gravity Field Refinement by Satellite to Satellite Doppler Tracking 
by Charles R. Schwarz 

December, 1970 

148 Inner Adjustment Constraints with Emphasis on Range Observations 
by Georges Blaha 

January, 1971 

150 Investigations of Critical Configurations for Fundamental Range Networks 
by Georges Blaha 
March, 1971 

177 Improvements of a Geodetic Triangulation through Control-Points 
Established by Means of Satellite or Precision Traversing 
by Narendra K. Saxena 
i June, 1972 

184 Coordinate Transformation by Minimizing Correlations Between Parameters 
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Preprocessing Optical Satellite Observational Data” 

3rd Meeting of the Western European Satellite Subcommission, IAG, Venice, 
Italy, May 1967. 

’’Global Satellite Triangulation and Trilateration” 
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’’Investigations in Connection with the Geometric Analysis of Geodetic Satellite 
Data” 

GEOS Program Review Meeting, Washington, D. C. , Dee. 1967. 
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50th Annual Meeting of the AGU, Washington, D. C. , April 1969. 
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* 

DEFENSE MAPPING AGENCY 

TOPOGRAPHIC CENTER 
WASHINGTON. O.C. 20315 


29 APR 1975 


Mr. Alfred Leick 
De>partment of Geodetic Science 
Ohio State University 
1958 Neil Avenue 
Columbus, Ohio 43210 

Dear Mr. Leick: 

Reference is made to your letter to Dr. Kenneth I. Daugherty, dated 
21 October 1974, requesting coordinates of Doppler stations in South 
America and your telephone conversation with Mr. John Love of this 
Center on 21 April 1975. 

The attached data partially fulfill your request. Enclosed for your 
retention are copies (front and back) of Geodetic Summary cards con- 
taining South American datum (SAD) positional data and Doppler Receiver 
Geodetic Summary sheets containing satellite derived positions for 
the following stations: 

STATION NO. ■ . STATION LOCATION 

General Conesa, Rio Negro, Argentina 
Villa Dolores, Argentina 
Tierra Del Fuego, Cerro Sombrero, Chile 
Frutillar Alto, Chile 
(Pre- and Post- Earthquake Values) 

Punta Arenas, Chile 
Santiago, Chile 
Arica, Chile 
La Paz, Bolivia 
Quito, Ecuador 
Coromoto, Venezuela 
El Callao/Tumeremo, Venezuela 

*1963 Provisional South Chile datum, not yet related to SAD. 

Authorization from the governments concerned has been obtained to 
release the coordinates of the Doppler stations for your use at Ohio 
State University; publication of the coordinates would require further 
authorization. 

jags bun* mqj 


30009 (T-009) 

30010 (T-010) 
*30011 

30012 

*30013 

30022 

30023 

30120 

30121 
30196 
30209 
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DMATC-G(52321) 

Ohio State University 


29 APR 1975 


Positional data for Station No. 30122, Asuncion, Paraguay, are not 
furnished because authorization for its release has not been obtained. 
Upon receipt of authorization from the Government of Paraguay, data for 
Station No. 30122 will be forwarded. 



Chief, Department of Geodesy 

Enclosure 
As stated 
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DOES MEAN SEA LEVEL SLOPE UP OR DOWN TOWARD NORTH? 


Comments on the article of the same title by Irene Fischer [1975] 


by 

M. G. Arur and Ivan I. Mueller 
Department of Geodetic Science 
The Ohio State University, Columbus 


1. Introduction 

Geodesists and oceanographers have disagreed on the direction and 
magnitude of the North -South sea level slopes along the East and West 
Coasts of the United States. 

There was some room to doubt the validity of the comparisons between 
the results of the geodesists and the oceanographers since they use different 
methods and different reference surfaces for the determination of these 
slopes [Fischer, 1975]. 

An attempt has now been made to compare the results of the ocean- 
ographers and geodesists by reducing them to the same terms* 


2. Method of Calculation 


The results of both geodetic and oceanic leveling have been reduced to 
the following compatible quantities for comparison at several stations along 
the two coasts (see Fig. 1): 

(i) Geopotential differences between the mean sea level and die deep sea 
isobaric surface used as a reference surface in oceanic leveling. 

(ii) Orthometric heights between the mean sea level and the, same deep sea 


isobaric surface. 

PRECEDING PAGE BLANK NOT FILMED 
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Values at the various stations for the anomalous geopotentials converted to 
dynamic heights of the mean sea level with respect to the standard ocean surface 
(0-decibar) of reference have been taken from the graphs of Sturges [1974]. Values 
at these stations for the orthometric height differences between the mean sea level 
and the reference geopotential surface used in geodetic leveling have been taken 
from the graphs of Balazs [1973]. 

Computations have been carried out for 21 stations along the TJ.S. East 
Coast and for 8 stations along the West Coast. The following assumptions 
have been made; 

(i) The deep sea isobaric surface used as a reference in oceanic 
leveling is an equipotential surface. 

(ii) Oceanic and geodetic leveling is in perfect agreement at Neah 
Bay on the West Coast and Port Maine on the East Coast, both 
having been used as references (origins) in geodetic leveling 
along the West and East Coasts respectively. 

(iii) The gravity field of the earth is well described by the normal 
gravity field in the areas under investigation and gravity varies 
linearly with height/depth up to 2 km . 

None of these assumptions will effect adversely the conclusions to be drawn. 

« 2.1 Calculation of Geopotential Differences 

2.11 Oceanic Leveling 


In accordance with the notation in Fig. 1, the difference in the geopo- 
tential between the deep isobaric surface and the mean sea level as deter- 
mined from oceanic leveling at an arbitrary station i , is 


AW a = AW + dW A * g X H + h X 1000 
0. s 0. m. s. D. 

l 1111 


( 1 ) 


where 


AW is the standard geopotential difference between the deep sea iso- 
s 


: * \ 
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baric surface and the standard ocean (0-db) surface. The value for 
this is 9704. 032 m 2 s -2 for 1000-db surface used as a reference on 
the West Coast [Montgomery, 1973]; and 19364.2 m 2 s' 2 for 
2000-db surface used as a reference on the East Coast [Bjerknes 
and Sands trom, 1910]. 

dWo i is the anomalous geopotential at station i, as determined by ocean- 
ographers . 

» Dl is the dynamic height at station i from the graphs of Sturges [1974] 
and is to be interpreted as per equation (1) above, 1000 being the 
constant value of gravity (in gals)used by Sturges to convert his original 
dW(j. potential anomalies into metric units. 

H S j is the orthometric height difference at station i corresponding to the 
standard geopotential difference AW S . 

g m . is the average normal gravity between the mean sea level and the deep 
sea isobaric surface at station i (in gals) based on the Geodetic Reference 
System. 1967. It is computed from 

g m . = 978.03185 (1 + 0. 005278895 sin^pj + 

H 0i 

+ 0.000023462 sin 4 <pi “ 0.000000227 7 )cms* 2 ( 2 ) 

2 

where <p j is the latitude and Ho^ is the orthometric height between the 
deep sea isobaric surface and the mean sea level (in meters). The constant 
0.0000002277 (x 978) is the normal vertical gradient of gravity in water. 

With the above notation AWqj can also be computed from 

AW 0i = g mi ^ H0i- (3) 

The quantities g m . and Ho i are to be determined iteratively until equations 

(1), (2) and (3) are mutually satisfied. 
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2.12 Geodetic Leveling 


The geopotential difference at the station i between the deep sea Isobaric 
surface and the mean sea level may also be computed from the results of 
geodetic leveling as follows: 

AW 0 - AW - dW (4) 

i g i 


where 

AW is the geopotential difference between the deep sea isobaric surface 
and the reference geopotential surface along which the geodetic 
’ leveling is assumed to take place, and computed at the reference 
station (origin of leveling) in accordance with section 2. 11. 

dW is the difference of geopotential between the reference geopotential 

g i 

surface of geodetic leveling and the mean sea level at the computation 
station, or 

dW s h x g , (5) 

CCS. 

B i B i i 

where hg. is the orthometric height at the computation station 
between the reference geopotential surface and the mean sea level, 
obtained from the graphs of Balazs [1973], and g S j is the average 
gravity along hg. obtained approximately by inserting Ho. - 0 
in equation (2). 

2.2 Results 

Figures 2 and 3 show the orthometric heights Hg. , i.e., the position 
of the mean sea level with respect to the deep isobaric surface as determined 
through oceanographic and geodetic leveling. The heights from oceanographic 
leveling were computed in accordance with section 2, 11, i.e., through an 
iterative, procedure to simultaneously fulfill equations (l), (2) and (3). The 
heights from geodetic leveling were computed from 
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H i - ’ 


( 6 ) 


H, 


where Hj was also computed iteratively as Ho i from oceanographic 
leveling, except using the potential difference AW instead of A Wq. . 

The differences between these two types of heights (oceanographic 
minus geodetic) near the mean sea level in terms of geopotential 
are shown in Figures 4 and 5.. 


3. Conclusions and Comments on the Paper by Fischer f 19751 

(i) The results of oceanic and geodetic leveling are comparable. 

In terms of magnitude, the discrepancies as pointed out by oceanographers, 
between oceanic and geodetic leveling, unfortunately do exist. 

(ii) If the deep sea isobaric surface is taken as the equipotential 
surface of reference, the results of both oceanic and geodetic leveling 
indicate that the ocean is sloping down from South to North along both the 
U.S. East and West Coasts. 

(iii) The discrepancies are greater along the West Coast where the 
deep sea isobaric surface of reference is 1000-db, as compared to the 
East Coast where the reference surface is 2000-db. The geopotential 
discrepancy is predominantly negative and increases with the distance from 
the reference station. 

(iv) Factors of some importance in the above comparisons are the 

lack of actual gravity data which made the use of normal gravity a necessity and 
the choice of deep sea isobaric surface. Both are likely to account for a very small 
part of the discrepancy in magnitude. The discrepancy about the direction 
of the slope seems to have been resolved. 

(v) Some of the above conclusions are unfortunately at variance with 
the findings of Fischer [1975]. The following comments are offered. 

In preparing her graphs (Figures 2, 3 and 4 of Fischer’s article), she 
has computed normal orthometric corrections (Ah) in ten-degree meridional 



sections from Helmert’s formula, multiplied by scale factors (k) so chosen 
that the sum of the corrections from the pole would yield an equatorial bulge 
of 2m, 1.6m and 1.8m, respectively, to match Sturges's profiles. According 
to Helmert’s formula the equatorial bulge for 100 m separation at the pole is 
53 cm, thus equipotential surfaces about 300 m apart at the pole will have a 
bulge of 1. 6m at the equator, which is the equivalent of a k factor of 3.019. 

Let us take the case of the West Coast and study the implication of consid- 
ering 1.6m as the equatorial bulge as represented in Figure 3 in Fischer’s 
article. This bulge and the use of the corresponding k factor (3.019) implies 
that she was in fact dealing with two equipotential surfaces which have an ortho- 
metric separation of about 301. 9m at the poles and 303.5 m at the equator. The 
upper surface is represented by the curve G in her graph and is thus referenced 
to the lower surface (let us call it 0), taken as a straight line about 300m below 
the ocean surface. Such a surface has no real meaning in the oceanographic or 
geodetic sense. It is too low to be considered as the theoretical 0-db surface 
(to which Sturge’s profiles are' referenced) and too high to be considered as the 
1000-db surface (which is about 990m below the ocean surface). Since Fischer’s 
G profile in fact refers to this 0 surface and Sturges’s profile to the 0-db surface, 
the convergence of these two surfaces needs to be considered, i.e. , when com- 
paring the two profiles the convergence needs to be computed and added to Sturges's 
values to refer them also to the 0 surface. 

Alternately, Sturges’s values can be plotted with respect to a correctly 
depicted theoretical 0-db surface which would show a curvature with respect 
to 0, and in fact would be practically parallel to the G surface, being within a 
few meters from it. 

If the plotting would have been done correctly as indicated above, it would 
have been obvious that the magnitude of the discrepancies between the ocean- 
ographic and geodetic leveling persists, and that the results of oceanographic 
and geodetic leveling have been correctly compared in the past. 

The main point thus is that geodetic leveling is done with respect to an 


G 


identifiable reference geopotential surface passing through the origin of 
leveling. Oceanic leveling provides dynamic heights with respect to the 
theoretical 0-db surface. These two equipotential surfaces, being within 
a few meters from each other, are practically parallel and they have a 
varying separation with respect to the deep ocean equipotential surface 
(such as the 1000-db surface) used as a reference in oceanic leveling. This 
varying separation is taken into account in the process of pceanic leveling. 
Thus for practical comparisons, geodetic and oceanic leveling provide values 
with respect to two parallel surfaces separated by the dynamic height of mean 
sea level at the reference station (origin) for geodetic leveling. 
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. 1 . 


Statement of Work 


Perform an error analysis, based on assumed sets of satellite borne 
transmitting equipment and ground receivers, to determine the optimum use of 
such systems in connection with the science that can be obtained from CLOGEOS 
measurements. 

2. Available Computer Programs 

To achieve the goals vide para. 1, main computer programs available at 
The Ohio State University are briefly described below. 

2. 1 Goddard Trajectory Determining System (GTDS) 

The Program GTDS, acquired recently at OSU from Goddard Space Flight 
Center, is extremely versatile and has numerous operating modes and capabilities 
[COSMIC, 74]. 

Space craft dynamics used in the program includes gravitational acceler- 
ation for Sun, Earth (up to 15 x 15 non-spherical field), Moon (up to 4X 4 non- 
spherical field) and all the planets, drag acceleration and solar radiation model 
with shadow effect and variations depending on distance from sun [Wagner, et. al. , 
1972]. 

In this project, this program was mainly used for ephemeris generation. 

2.2 The Ohio State University Geometric and Orbital (Adjustment) 

Program (OSUGOP) 

The program OSUGOP performs basically as an adjustment program, using 
optical or range observations, in geometric /orbital mode [Reilly, et. al. , 1972], 
One important feature of the program is its capability to apply and obtain different 
constrained solutions, including 'inner' constraints LBlaha, 1971]. 

As the current error analysis was to be based on the ’estimable' quantities 
e.g, , chords and angles between the locations of ground stations, a new subroutine 
'CHECK' has been added to the program. 
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2.3 Short Arc Geodetic Adjustment Program (SAGA) 

, The latest version of this program, as developed by Duane Brown Associates, 
was obtained at OSU from Air Force Cambridge Research Laboratories [Brown 
and Trotter, 1973]. 

The program employs a power series solution using partitioned regression 
technique. Inner constraints capability has also been added in the present program. 

2.4 Range Generation-Geometric - with 7 Records in Mode B 
Program (RGGR7-B) 

The computer program RGGR7-B generates ranges from the short arcs on 
tape as generated by GTDS and given station locations. Time interval between the 
ranges, cut off angle or maximum zenith angle as well as the station locations and 
the model in which they observe the satellite per pass can be specified outside the 
program. The program has the capability to superimpose white noise of any 
standard deviation on the ranges. The ranges are written on tape in the GEOS 
range format as required for the input to the computer program OSUGOP. 

3. Data Generation 

3. 1 Satellite Orbits 

During the orientation meeting for CLOGEOS at Huntsville, Alabama on 
February 6, 1975, it was agreed upon that the possible approximate altitudes for 
investigation are 350 km, 600km and 1000 km. However, during the period under 
report some alterations were made in the designated orbits. 

The orbit generation was carried out in two steps - first a long arc for 
126 hours was generated in inertial system and then taking the suitable orbit 
points several short arcs in the body fixed system were generated over the area 
under investigation. 

The details about orbit delineation and the computer expenditure involved 
are given in Tables 3 . 1 - 1 to 3 . 1 -5 . 

3.2 Range Generation 

The computer program RGGR7-B and the short arcs (vide para . 3.1) were 
then used to generate ranges in geometric mode with Gaussian standard deviation 
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of 10 cm. During the range generation the density of orbit points was suitably 
altered between lower, middle and upper' orbit to keep the number of ranges the 
same for each orbit. 

During the range generation a new variable was introduced as under (Fig. 1): 

Case A - When all the observing stations were located on the ellipsoid 
(h=0). 

Case B - When some of the observation stations were raised from the 
ellipsoid to a maximum height of 100 meters . 

Case C - Same as Case B except maximum height was made up to 
1000 meters. 

The station layout, orbital parameters, span of observations, general dis- 
tribution, etc. are given in Fig. 2 and Fig. 3. For details of data generated and 
the computer expenses see Table 3.2-1. 

The ranges were generated on three tapes RNGE01, RNGE02, and CSTPOl 
and the contents of files on these tapes are given in Tables 3.2-2, 3. 2-3 and 3.2-4. 

Tape Format 

DCB=(RECFM=VBS, LRECL=40, BLKS IZE=12004) 

Record Format 

NN, IYMD, IH, IM, IS, RR, SIGR 

where 

NN = Station identification number 
IYMD = Year, month and date of observation in packed format 
IH = Hour 
IM = Minute 
IS = Seconds * 10 4 
RR= Range in meters 

SIGR = Standard deviation of range in meters. 
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4. 


Simulated Solutions 


The CLOGEOS error analysis was decided to be carried out in two main 
section viz. , geometric and short arc modes. During the period under report ' 
the investigations were made only in geometric mode. Thirty-seven simulated 
solutions were computed as detailed in Tables 4-1 to 4-3. 

5. Analysis and Conclusions (Preliminary) 

In an error analysis of the type under consideration it is more realistic 
to analyze the estimable quantities. As the station coordinates fall under non- 
estimable quantities, error analysis for them would have been significantly 
dependent ‘on the origin and the error propagating outwardly from it. Distances 
and angles are the only estimable quantities in the present error analysis and 
are defined as under: 

chord R tJ -x/+ (y, -y s f+ (z\ - Zj f 


, (Xi-x J )(x k -x J )+(y 1 -y 1 )(y k -y J >f(z I -Zj)(z k -z, ) 

angle Q! t Jlc = cos 1 ===== t === ■■■- • === 

/ ( x i -Xj f+(Yi -Yi f +( z i -Zj f / (x k -x. f-Hx-yj f+( z k -z, f 


(5.2) 


where x t , y t , z } are the rectangular coordinates of the i points. The current 
investigation deals with the distance analysis and the results have been broadly 
grouped as under. 

5. 1 Effect of Orbit And /Or Station Separation 

Six solutions were run for 500 events each for cases A, B and C and the 
standard deviations for the best case o r (where r^ denotes the distance between 

3 2 

station numbers i and j) and for the worst case a r were plotted (Figures 4 and 5). 

19 

The a - s show significant improvement either introducing station 
r U 

separation as in Case C or by introducing mixed satellite altitudes through 
eccentric or multiple arcs . . 
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5.2 Effect of Orbital Height 


Figures 6 (for a ) and 7 (for a ) show that the height of the orbit 

r i2 r i 9 

cannot improve solutions when observing stations have a ’near' critical configu- 
ration, 

5.3 Effect of Number of Events 

Figures 8 (for a ) and 9 (for ct ) show that the number of events 
r ia r i s 

included in any solution had a significant influence. However, in any given 
station configuration, the less critical is the configuration the smaller the number 
of events needed to obtain a specified accuracy. 

6. Personnel 

Ivan I . Mueller, Project Supervisor, part time. 

Muneendra Kumar, Graduate Research Associate, part time. 

Boudewijn H.W. Van Gelder, Graduate Research Associate, part time. 
Michelle Neff, Administrative Assistant, part time. 

7. Travel 

1. Project meetings at MSFC, Huntsville, Alabama. 

December 17, 1974 (Mueller and Van Gelder). 

February 6, 1975 (Mueller, Van Gelder and Kumar). 

See material distributed in Appendix A. 

April 17, 1975 (Mueller) 

Material presented at this meeting is identical 
to the one in this report. 

2. January 30 - Feb. 2, 1975, New York City, (Mueller) 

To attend the AAAS meeting. 

3. April 2-5, 1975, Siena, Italy, (Mueller) 

Mathematical Geodesy (Partial support). For a 
report on this meeting see Appendix B. 
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Table 3.1-5 


Short Arcs 


Satellite 

Height 

(km) 

No. of 
Passes 

Length of each 
Pass in Time 

Density of 
Satellite points 

Computer Expenses 
$ 

9 

■■ 


10/sec 

300.00 

392 

35 

8 min 

10/sec 

1050.00 

657 

22 

10 min 

10/sec 

770.00 

1007 

30 

12 min 

1 0/sec 

1200.00 


3320.00 
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Table 3.2-1 


Satellite 

Height 

(km) 

No. of 
Passes 

Case 

Type 

Maximum 
Data Points 
Generated 

Computer Expenditure 
$ 



A 

3000 

80.00 

9 

30 
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1500 

60.00 
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1500 

60. 00 



. 
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5000 
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392 

35 
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500 

30.00 
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500 
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A 

500 
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657 

22 
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500 

30.00 

. 
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500 

30.00 



A 

5000 

180.00 

1007 

30 
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5000 

180.00 
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500 

30. 00 

I 

I 
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$920.00 
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Table 4-1 


Satellite Orbits Varied 


Case 

No. /Type 

No. of events 

Total 

No. of 

Computer Expenses 

Type 

of orbit 

per orbit 

Events 

Simulated 

$ 


used 

• 



Solutions 



. 

l(low-low) 

500 

500 

3 


A,B 

l(lower) 

500 

500 

3 


and C 

l(middle) 

500 

500 

3 

18x50.00= $900.00 


1 (upper) 

500 

500 

3 



2 ( lower + 

250 

500 

3 



middle ) 

' 





3(lower+ 

133 

500 

3 


■ 

middle + 
upper) 















Table 4-2 


No. of Events Varied 


Case 

Type 

Type of 
Orbit Used 


No. 

of events included in the 
Simulated Solution 

Total* 
number of 
Solutions 

Computer 
Expens es 
$ 

■ . 

low - low 

50 

500 


1500 

1 

^^2 


3 

120. 00 


lower 

50 

500 

1000 

■ 

2500 

■ 

5000 

4 

400. 00 

A 

middle 

50 

500 

1000 

■ 

2500 

- 

5000 

4 

400. 00 


upper 

50 

500 

1000 


2500 

- 

5000 

• 

4 

400. 00 

B&C 

low -low 

50 

500 

■ 

— 

1500 

- 

„ 


4 

80.00 


$1400.00 


*Solution for 500 events excluded as already counted in previous Table. 
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(ii) Ice-flae Monitoring 


(iii). Iceberg Monitoring 


s (iv? Glacier Monitoring 
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6V those Points i.e. local 

COOR3MNBTE SVSTEM 

SYSTEM RRNK ^DEFICIENT “J CQ 

USING THE <3ENE ^RUZEO INVERSE 

90LOT10N £\K)NER CON^»TRB\KiT ^ 
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: NOW ESTIMABLE 

QOBNTmES 


THE OHIO STATE UNIVERSITY 
DEPARTMENT OF GEODETIC SCIENCE 
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Appendix B - 1 


Trieste, February 18, 1975 

6th SYMPOSIUM ON MATHEMATI CAL GEODESY 
(3rd Hotine Symposium) 

Siena, April 2-5, 1975 


Dear Colleague, 

Following to my letter of January 8, I am sending the 
Second (and last) Circular Letter concerning the 6th Symposium 
on Mathematical Geodesy (3rd Hotine Symposium) organized by 
the Italian Geodetic Commission under the aegis of S.S.G. No. 

4*31 "Mathematical Techniques in Physical Geodesy" . 

The opening of the Symposium will take place at 9.30 
A. M. of April 2. 

The participation to the Symposium is not bound to any 
formality or payment of a registration fee . You are only re- 
quired to fill in the enclosed reservation form, and send it 
directly to the J.zienda Autonoma Turismo, Siena, as indicated. 

I recommend Hotels Minerva, Continental, Toscana, Chiusarelli 
and Pension Ravizza, that are close to the University in which 
the Symposium will take place (Via Banchi di Sotto 57, 1st floor). 

The dead line for the hotel reservation is March 10; an 
earlier application will be highly appreciated. 

In line with the former Hotine Symposia, no preliminary 
presentation of papers is required; however I will appreciate 
receiving the title of your contribution. The main accent will 
be on oral presentations and discussions on unusual aspects of 
Geodesy. The publication of the Proceedings of the Symposium is 
foreseen. 

Looking forward to seeing you in Siena, and with my 
heartiest regards and wishes . 

r 

Yours very sincerely, 

/ 

prof .Antonio Marussi 

flh • 

PRESS'D LMSrirU'IO DI GEODESIA L GEOFISICA DELL' II M1VE.RSITA DEGL1 STUD! DI THUSTc. 
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March 17, 1975 


Professor Antonio Marussi 

Direttore deli'Isiituto di Geodes ia e Geofisica 

Bell'Universita di Trieste 

Via deU'Universita 7 

34100 Trieste 

Italy 

Dear Professor Marussi: 

This is just a short note to let you know that I will be present at 
the Kotina Symposium in Siena. I do not plan to present a formal 
paper, but I would like to have some private discussions on the 
unusual aspects at the Close Grid Geodynamics Satellite Measure- 
ment System (CLOGEOS). This proposed system consists of a 
number of satellite-borne lasers and closely spaced reflectors 
on the ground for die primary purpose of monitoring fault motions, 
etc. 


I am looking forward to seeing’ you again, 
to you and Mrs. Marussi 


With my best wishes 


Sincerely yours. 
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6th SYMPOSIUM ON MATHEMATICAL GEODESY 
(3rd Hotine Symposium) 

Siena, April 2-5, 1975 

PROGRAM 

Wednesday - April 2, 1975, 15-17 
Chairman: I. Mueller, Columbus/U.S. A. 


1. Cartan and the holonomity problem I 

A. Marussi, University of Triests, Trieste/Italy 


2. Cartan and the holonomity problem II 

N. Grossman, University of California at Los Angeles 
(UCLA), Los Angeles /U.S. A. 


3. Cartan and the holonomity problem III 

F. Bocchio, University of Triests, Trieste/Italy 

4. Cartan and the holonomity problem IV 

E. Grafarend, University of Bonn, Bonn/W- Germany 

5 , Cartan and the holonomity problem V 

J.G. Leclerc, University of Stockholm, Stockholm /Sweden 
and Quebec /Canada 
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Thursday, April 3, 1975, 9-l2h 

.... ,0 

Chairman: T. Krarup, Copenhagen/Denmark 

1. Utilization del documents cartographiques existants (anomalies de Bouguer, 
cartes d ’altitude) pour une definition precise du potentiel dans l'espace 
exterieur au geoide vrai 

H.M. Dufour, IGN Paris /France 

2. On the determination and application of gravity gradients in geodetic systems 
E. Groten, Technical University of Darmstadt, Darmstadt/Germany 

3. Approximation of certain solutions of the exterior oblique derivative problem 
for the laplace equation 

K.J. Witsch, University of Bonn, Bonn/W -Germany 

4. Boundary problems for the sphere 

E . Ecker, Technical University of Berlin, Berlin/W -Germany 

5. Analytical continuation of a function from the length’s surface upwards 
M. Pick, Academy of Science, Prag/CSSR 

Thursday, April 3, 1975, 14-l7h 
Chairman: E. Grafarend, Bonn/Germany 

1. Reflexive predictions 

A. Bjerhammar, University of Stockholm, Stockholm /Sweden 

2. Determination of datum-shift parameters using least-squares collocation 
and 

A mass density covariance function consistent with the covariance functions 
of the anomalous potential 

C. Tscherning, Geodetic Institute Copenhagen/Denmark 
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3. A spherical harmonic expansion of the isostatic reduction potential 
G. Lachapelle, Geodetic Survey of Canada, Ottawa/Canada 

4. Least squares collocation for large systems 

K.P. Schwarz, Technical University at Graz, Graz/Austria 

Friday, April 4, 1975, 9-l2h 

Chairman: A. Bjerhammar, Stockholm/Sweden 

1. Unusual aspects at the close grid geodynamics satellite measurement system 
(CLOGEOS) 

I. Mueller, The Ohio State University, Columbus /U.S. A. 

2. Free adjustment of a torsion balance net 

G. Hein, Technical University of Darmstadt, Darmstadt/Germany 

3. The nature of space near the earth 

N. Grossman, University of California at Los Angeles (UCLA), 

Los Angeles, U.S. A. 

4'. A general method for the computation of minimax -errors 

G. HeindlandF. Reinhart, Technical University of Munich, Munich/Germany 

Friday, April 4, 1975, 14-17h 

Chairman: H.M. Dufour, Paris/France 

1. 3-d Mapping and mapping of the gravitational field 
A. Marussi, University of Trieste, Trieste/Italy 

2, On a general potential- invariants representation of the geopotential field and 

applications, ‘ 

F. Bocchio, University of Trieste/Italy «. ' • 
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PREFACE 


This project is under the supervision of Professor Ivan I. Mueller, Department 
of Geodetic Science, The Ohio State University and is under the technical direction of 
Dr, Nicholas C. Costes, Code ES31, MSFC, Huntsville, Alabama. The contract is 
issued by Procurement Office, MSFC and is administered through ONRRR, Columbus, 
Ohio. 


1 . 


Statement of Work 


Perform an error analysis based on assumed set of satellite borne transmitting 
equipment and ground receivers to determine the optimum use of such systems in 
connection with the science that can be obtained from CLOGEOS measurements . 

2. Data Generation 

2. 1 Range Generation 

The computer program RGGR 7-B (as used for geometric mode) was suitably 
modified to generate ranges for analysis in short arc mode. Using the short arcs 
(see paragraph 3.1 of First Quarterly Status Report), ranges were generated with 
a Gaussian standard deviation of 10cm. One set of errorless ranges was also 
generated. 

For details of data generated and the computer expenses see Tables 2.1-1, 

2. 1-2 and 2.1-3. 


3. Simulated Solutions 

During the reporting period the investigations v/ere mainly made in short arc 
mode. Thirty-one simulated solutions were computed as detailed in Table 3-2. 

In addition, twelve simulated solutions (Table 3-3) were also computed in 
geometric mode. Computer expenses for simulated solutions are given in Table 
3-1. 

4. Analysis and Conclusions (Preliminary) 

4.1 Geometric Mode 

Reference to equations 5. 1 and 5.2 of the First Quarterly Status Report, 
some further graphical analysis was performed to study the effect 


- 2 - 


of orbit, station separation and number of events. Only the "typical” angles a ljk 
and the distances - r^ were used in this analysis. 

The above graphical analysis gave the same conclusion for the geometric mode 
as reported in the First Quarterly Status Report. However, the details will be included 
in the final report. 

4. 2 Short Arc Mode 

4. 2. 1 Effect of Fundamental (LAGEOS) Stations Per Pass 

Figure 4. 2. 1-1 shows the effect of inclusion of fundamental (LAGEOS) stations 
in short arcs. The reduction of LAGEOS stations to two per pass adversely effects 
the solution and for the directional stability, each pass must have at least two funda- 
mental stations. 

4.2.2 Effect of Grid (CLOGEOS) Stations Per Pass 

Figure 4. 2. 2-1 and Tables 4.2.2-land 4. 2. 2-2 show the effect on recovery of 
relative position of grid (CLOGEOS) stations as the number of such stations varies in 
any satellite pass. 

4.2.3 Effect of Orbital Height 

Table 4 . 2. 3-1 and Figures 4.2. 3-1 and 4.2. 3-2 show the effect of varying orbital 
height of the satellite. Even though the residuals in coordinates for the lower orbit 
(Table 4.2. 3-1) are larger compared to upper orbit case, the overall recovery (Figures 
4. 2. 3-1 and 4.2. 3—2) is quite compatible in both cases. 

4.2.4 Effect of Number of Events 

Table 4 . 2. 4-1 and Figures 4 .2.4-1 and 4. 2.4-2 show the effect of different 
number of events in a simulated solution. The results show that variations are not 
significant. 
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4.2.5 


Effect of Observational Mode 


Figure 4. 2. 5-1 shows the effect of simultaneous (grouped) versus the sequential 
observational mode in the short arc mode. The recovery is more or less the same in 
each case. 

4. 3 Geometric Mode vs . Short Arc Mode 

Figures 4. 3-1 and 4.3-2 show the comparative recovery in residuals in geometric 
(with 9 stations) and short arc mode (with 12 stations) for cases A and C. 

In case A, the residuals for geometric mode are in meters, while all other residuals 
are in centimeters. 

Figure 4.3-3 and Table 4. 3-4 show the comparative recovery in residuals where one/ 
three LAGEOS stations have been added in the geometric mode also. 

5. Personnel 

Ivan I. Mueller, Project Supervisor, part time 
Muneendra Kumar, Graduate Research Associate, part time 
Boudewijn H. W. VanGelder, Graduate Research Associate, part time 
Michelle A. Neff, Administrative Assistant, part time 

6. Travel 

1. June 15 - 20, 1975, Washington, D* C. (Kumar) 

Annual Meeting of American Geophysical Union 
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Table 2.1-1 


Long Arcs 


Satell ite 

Oribt 

Time of Data 

Height 

Type 

Generation 

(km) 


(hrs) 

392 

Polar 

126 

657 

-do- 

126 

1007 

-do- 

126 


Short Arcs 


Satellite 

Height 

(km) 


No. of 
Passes 


Table 2.2-2 


Length of each Density of 
Pass in Time Satellite points 



8 min 


10 min 


12 min 


l/sec 


1/sec 


l/sec 


775.00 


900.00 


2,215.00 
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Tor Nni - 



€>R»D 

STATIONS XJRYS 
PER 


EVENTS , RESULTS 

Hw.euItotrl totrl ^ RESI-DUBLS INI CM 


P. PASS EVENTS OBSSRVfTfT. 



5 UPPER o S 3 9 3 9 S 2.39 1 S &8 \i 10 


8 Lower 



5 7 lo 254 3048 3 


2 UPPER •. « 


13 Lower •• 



3 9 15 983 5796 6 


5 7 29 758 9^96 







































































Thbue ^.2.^ - 1 


s R ESULTS 

CM ]RES\3>0<=iLS IN CM 


S.VENTS OCSEKvqT. 



5 upper O S 3 9 3 9 8 239 2.848 12 10 


« IS 953 5795 L 




5 16 » 969 1162s •> •• 


» 7 10 259 3048 3 .. 


.. .. 29 758 9096 















































































GEOMETRIC MODE vs. SHORT FIRC MODE, 


# ORBIT 



FUNX^J GRID EVENTS 

neNTfiL STFvnoNS -mwc. — — — i , /\r 

sTfmows per pr^s T ° TRL 

nor£ PER PASS PHAA P-PflSS EVEOTS 


1 


RESULTS 


RESl^UfiLS IN CM sFKtftt CT 

cn Cf1 


S.A. 2. upper o & 3 9 3 9 15 ^83 \Sl^ 6 ho 


a 


G. 1-9 


<X 

■E-7 


'* W 


■i 

■iuhu 





■ 



507 hsto h * 




■ 


" 6i8o •• 


•• ^83 579^ 




2 3 

^ 2 

























































































STRTIONS 


GRID STRTIONS 

1005 <P = 38° O 1 X = 2^0° o* 

OTHER. & GR\T> QTfVtIOHS RRE. 6«OoP&» RR00N3> \005 RT 5* 


FUNDRMENTRL stations 

IOIO SRN 35 IESO M> S 33 ° o' X S 2 ^ 3 * 0 * 

lot I ©ERR LAKE = ^2° o' s 2.^8° 30* 

1052. Quincy « ^0* o* s 239° o 1 


interval a 




£FFjrQTT OF NUMBER OF G>R\3? STRTIONS PER 


U PPER ORfc\T \ \OOCj KM 

» • / 

250 EVENTS \ 

' ..... » 

CASE. A 

1 * 

3 funx>bh. stftt ions'* PER /prss 

\ # 

•;;; L\ GR iX> \ ST./ PER PRSS 
~ 6 G\-UD V. PER PRSS 
Cj G R|»B ST. PER PRSS 


» \ 





Fie - \ 


EFFECT OF OR&ITRL HEIGHT 


2.50 EVENTS (3000 OBS.) 
HEIGHT FRCTOR O 


LOWER ORBIT 392. KM 


UPPER ORBIT VOO-7 KM 



EFFECT OF ORB1THL HEIGHT 

500 EVENTS (6000 OBS.) 

Height frctor O 

LOVOER ORBIT 3^2. KM 

UPPER ORBIT lOOT KM 



F\e kl.3 -2. 







20 T f'SJCM 


EFFECT OF OBSERVATIONAL MocdE 


F\G 4.2.5 - | 


UPPER ORG>\T \007 KM 

HEIGHT FRCTOR O 

250 EVENTS (3000 OBS.) 


GROUPED obsgrvrti oms 

5EQOENTIRL OSSERVRTVOHS 


1 z 1 *1 *• i 'll 



GEOMETRIC MCKDE VS. 6 


UPPER ORGrr Voocj km 
500 EVENTS 



GEotnETR\c Hoa>e 




SHORT PRO MODE 



NOn^ER OF RE3D0t=»-S Vrv* 


Hetight 

IRRO HCCDE \ ppjctor ^ 


FiS ^.3-1 


LOWER ORGIT 3 CJ 2 . Kh 
500 EVENTS 

< 9 > 





GEOtlETRYC. MODE 


r — j 


r 


GHORT PRC MODE 



GEOMETRIC MODE \/S. 


UPPER ORBIT \OOCj Kh 
500 EVENTS 

GEOMETRIC MODE (^) 
"1 SHORT RRC MODE 



HEIGHT 

SHORT RRC MODE ; p RCT0R 1000 


P\e A3 -2 


LOWER ORBIT 3(^2 KM 
500 EVENTS 

geometric mode (9) 

SHORT RRC MODE , CASE R 






GEOMETOC HoDE 


SHORT F1RC MO30E1 


F\g ^. 3-3 


UPPER ORB\T 


too^ 


500 BVE.i\rrs 


GEOMETRIC 
Mo DE 

C) 6TRTIONS 


r— ^ I 


GEOMETRIC 

MODE 

12 STRTIONS 


SHORT RRC 
... MODE 
! 12 STATIONS 


Height ffictor o 


Height factor looo 


$ 

u 

u 

\A 

Vy 

.... M 
I 0* 

‘••f / A 

\/A 

r? / •' 


I I I * t I o iu *- 

jc- w *» — .o (n •*, tn 

'«/> Vt iff" VI V\ VI 


</> 

J 

CE 

10 2 


in O' j- 

* ^ Icn 


i - ✓ 

i 

i • • 

r~: 

t : ■ 
f : 

. J 


-■ - * J — T~ 1 1 - 

i i ». i • < i i o r w Jt o> 

^ ^ * ? In ^ t/1 «/| tfl «/> (/I 

*n V) V\ in V\ ij\ ^ 1 1 


S'SZ 



